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ABSTRACT 

We present a simplified model to study the orbital evolution of a young hot Jupiter 
inside the magnetospheric cavity of a proto-planetary disk. The model takes into ac- 
count the disk locking of stellar spin as well as the tidal and magnetic interactions 
between the star and the planet. We focus on the orbital evolution starting from the 
orbit in the 2:1 resonance with the inner edge of the disk, followed by the inward and 
then outward orbital migration driven by the tidal and magnetic torques as well as the 
Roche-lobe overflow of the tidally inflated planet. The goal in this paper is to study 
how the orbital evolution inside the magnetospheric cavity depends on the cavity size, 
planet mass, and orbital eccentricity. In the present work, we only target the mass 
range from 0.7 to 2 Jupiter masses. In the case of the large cavity corresponding to 
the rotational period ^ 7 days, the planet of mass > 1 Jupiter mass with moderate 
initial eccentricities 0.3) can move to the region < 0.03 AU from its central star in 
10^ years, while the planet of mass < 1 Jupiter mass cannot. We estimate the critical 
eccentricity beyond which the planet of a given mass will overflow its Roche radius and 
finally lose all of its gas onto the star due to runaway mass loss. In the case of the small 
cavity corresponding to the rotational period ^ 3 days, all of the simulated planets lose 
all of their gas even in circular orbits. Our results for the orbital evolution of young hot 
Jupiters may have the potential to explain the absence of low-mass giant planets inside 
~ 0.03 AU from their dwarf stars revealed by transit surveys. 
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Introduction 



The number of discovered extrasolar planets has been incre asing; quickl:\{l| since the first extra- 
solar planet, 51 Pegasi b, was detected around a solar-type star (IMavor fc Quelodll995l ). 51 Pegasi 
b is the prototype of a class of exoplanets known as "hot Jupiters" because they are generally 
Jupiter-mass planets located ^0.1 AU from their parent stars and therefore their equilibrium tem- 
peratures are much higher than that of Jupiter. The high temperature of the protoplanetary disk at 
the orbital radii where they are now provides an environment ex tremely different from where g:iant 



plane ts are expected to form via either the co re-accretion model (IPollac 



2002) or the gravitational instability model (lGammiell200ll : IBoss 



c et al 



2004 



1996 



Kokubo fc Ida 



RafikQv|l2005|). One of the 



most commonly adopted solutions is that hot Jupiters have formed at larger orbital radii in the pro- 
topl anetary disk and then moved inward to their current locations via the disk-planet interactions 
(e.g. 



Lin et al. 



19961). 



The magnetic fields of Classical T Tauri stars (hereafter CTTS) are typically strong enough 
to truncate the inner regions of protoplanetary disks to the corotation radius, where the Keplerian 
angular velocity of the disk is the same as the stellar-spin angular velocity, and create an inner 
magnetospheric cavity. Besides, the migration due to the planet-disk interaction is ex pected to slow 
dramatically once a young hot Jupiter passes the inner disk edge ( Lin et al. 19961 ). The pile- up 
population of hot Jupiters at ^ 0.0 4 AU from their solar-type parent stars revealed fr om radial- 
velocity searches (iGaudi et al.ll2005l ) may be attributed to the existence of the cavity (ICarrI 120071 . 
and references therein). 



On one hand, an atte nipt has been niade bylSetiawan et al.1 (120081 ) to search for a hot Jupiter 
around a CTTS, although iHuelamo et all (120081 ) substituted an explanation of a cool stellar spot 
for an orbiting hot Jupiter. Future observations will provide more evidence for planets of early ages. 
On the other han d, the simulated o rbital evolution of giant planets as they migrate into the cavity 



was carried out by lRice et al.l (120081 ). Their results suggest that high-mass hot Jupiters would move 
closer to their central stars and even be destroyed because planet's entry into the magnetospheric 
cavity results in a growth in orbital eccentricity e. Although the eccentricity excitation is subject 
to the uncertain properties such as disk mass and viscosity, this model may provide a mechanism to 
pump up the eccentricity of a hot Jupiter inside the cavity and in turn affect the orbital evolution 
of the planet via the tidal interactions between the star and planet. In addition, the magnetic 
interactions between the star and planet in the magnetospheric cavity would be important as well, 
which m ay be an ups cale analogy to the magnetic interactions between Jupiter and the Galilean 
satellites (IZarkal 120071 . and references therein). 



The coupled ti dal evolution of the radius and o rbit of hot Jupiters in t he absence of a disk has 
been investigated (IGu et all 12003 : iMardlingI 120071 : llbgui fc Burrowsl l2009l ) . To obtain noticeable 



influence on i?^, the tidal heating is assumed to be deposited deep in the planet in these models 



^ See the Extrasolar Planets Encyclopaedia (http://exoplanet.eu') for detail. 
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fjGu et alJbood ). The purpose of iMardlingl (120071 ) and llbgui fc Burrowsl (120091 1 is to explain the 
large radii of some transiting hot Jupiters using tidal heating. On the other hand. lou et al. ( 2003 ) 
showed that for modest eccentricities, young hot Jupiters of radius ^ 2 Jupiter radii at < 0.04 AU 
swell beyond two Jupiter radii and their internal degeneracy is partially liftecil. Thereafter, their 
thermal equilibria become unstable and they undergo runaway inflation until their radii exceed 
the Roche radius. Then Roche lobe overflow ensues. Under the assumption that the overflowing 
gas immediately loses its orbital angular momentum to the planet and plunges into the central 
star, these mass-losing planets migrate outwards, such that their semi-major axes and Roche radii 
increase while their mass, eccentricity, and tidal dissipation rate decrease until the mass loss is 
quenched. Of course, including a d isk would change the orbital evolutions of Roche-lobe-overflowing 
planets, as has been calculated by iTrilling et al.l (Il998l ) who do not consider the tidal expansion of 
Rp and the magnetospheric cavity though. 

As more and more hot Jupiters located < 0.04 AU have been discovered by transit surveys 
for recent years, it is noteworthy that almost no giant planets of mass < 1 Jupiter mass (Mj) 
have been detected within ^ 0.03 AU from their solar-type parent stars (see Figure [U also see 



Udry fc SantosI 120071 . and references therein). This puzzle may be attributed to the subsequent 
evolutions of hot Jupiters afte r they enter the magne tospheric cavity of the protoplanetary disk. 
Since thermal eva poration (e.g. iDavis fc \y heat levi 12009 ) has not been undoubtedly identifled to be 
the major cause (IMurray-Clay et al.ll2009l . and references therein), we shall discuss the potential 
implication of our model to explain the observation in the end of the paper. 

For the above reasons, in this work we concentrate on the orbital evolution of a young hot 
Jupiter inside the magnetospheric cavity of a protoplanetary disk by taking into account the tidal 
and magnetic interactions between the star and planet. To restrict to the problem without sig- 
niflcantly involving the planet-disk interaction for simplicity, we consider the evolution after the 
hot Jupiter enters the 2:1 orbital resonance with the inner edge of a disk. The initial eccentricity, 
which is expected to be excited in the cavity before the planet migrates to the resonance location, 
is treated as one of the free parameters of the problem. The equations for the orbital and interior 
evolution of a hot Jupiter inside a magnetospheric cavity are described in §2. The results are 
presented in §3, and are summarized and discussed in §4. 



2. Model Description 

Figure [2] illustrates our toy model that shows a hot Jupiter orbiting a CTTS in a magnetospheric 
cavity of a protoplanetary disk. The magnetic dipole flelds of a protostar truncate the inner part 
of the protoplanetary disk at the disk radius Rmdiski which is the radius of the cavity. The initial 



^For comparison, the initial orbits and eccentricities of a young hot Jupiter in Ibgui Burrows (|2009|) are larger. 
Therefore the tidal heating in their model can proceed at a modest rate to attain the large Rp on the timescale of 
Gyrs. 
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orbit of the planet is set to lie at the 2:1 orbital resonance with the in ner edge of the disk, which 



allow s us to ignore the tidal interactions between the disk and the planet (ILin et al.lll996l : lRice et al 



20081 ) . The semi-major axis a of this resonant orbit is denoted as a2.i- However, when the inner 



edge of the disk extends inwards for any reasons in our model such that a > a2:i, the planet is 
artificially moved inwards as well to maintain the 2:1 orbital resonance with the inner edge of the 
disk. Moreover, we assume that there is no magnetic linkage between the disk and the planet as 
illustrated in Figure [2] and therefore neglect magnetic interactions between them. In other words, 
the planet's orbit evolves due to tidal and magnetic interactions with the central star, whereas 
the star interacts with both the planet and the disk. Therefore, this section consists of two parts, 
namely the magnetic and tidal effects, to describe the model in detail. §2.1 and §2.2 describe 
the star-disk and star-planet magnetic interactions, which leads to disk locking and contributes to 
the orbital evolution of the planet, respectively. §2.3 and §2.4 then describe the tidal interactions 
between the CTTS and the planet. The circularization of the planet's orbit is due to tides on both 
the planet and the star. The thermal inflation of the planet is due to tides on the planet excited by 
the star. On the other hand, the tides on the star raised by the planet affect the evolution of the 
planet's orbit and the stellar spin. Finally, we summarize the key formulae for our model in §2.5. 



2.1. Star-Disk Magnetic Interaction 



There have been a number of models for the star-disk interactions to account for the rotation of 



Lovelace et al 



CTT Ss as well as star and jet formation (e.g. IShu et al.lll994l : iMatt fc Pudritdl2005l : 
20081 ) . In addition to the closed field lines connecting from the protostar to the disk to allow for 
funnel accretion, t he open fie l d line s may diverge primarily from the inner edge of the disk (i.e. 
X-wind model, see IShu et al.Ml994l) or emanate larg!;ely from the star and the disk (i.e. stellar 
and disk wind model, see iMatt fc Pudrit d l2005l). Furthe rmore, the stellar field topology may 
be in the complex form of multipole (ILovelace et al.ll2008l ) or tilted dipole relative to the stellar 
spin ( Laine et al. 20081 ). In this work, we ignore these complications and adopt the simple non- 
tilted dipole configuration of stellar fields that thread through the protoplanetary disk as has been 
illustrated in Figure [2 

To calculate the mag netic torque Tmqg on the s tar due to the magnetic linkage to the disk, 
we adopt the approach by lArmitage fc Clarkd (119961 ) in which the magnetic torque acting on an 
annulus of the disk is due to the azimuthal twisted field lines B^h caused by t he difference between 
the stellar spin and the orbital motion of the annulus (ILivio fc Pringlelll992l ). Hence, the torques 
integrated from the disk inner edge at Rmdisk to infinity contribute to the net magnetic torque 
between the star and the disk: 

,2 

(1) 



T — _ 9P-3/2p-3/2 x 



where ji — B^R^ is the stellar dipole moment, B^ is the magnetic field strength at the stellar 
surface, i?* is the stellar radius, Rc — {GM^/Vtiy/'^ is the corotation radius, fi* is the stellar spin 
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r ate, and is the stellar mass. We employ the same scaling law as used in lArmitage fc Clarke 
(jl996l ) to specify the strength of stellar magnetic field: 

where the scaling factor Bq is a constant indicative of the surface magnetic strength for the stellar 
spin period of 4 days. 

The magnetospheric radius Rmdisk is determined by the balance between the stellar magnetic 
stress and the disk Reynolds stress (Starczewski et al. 2007, and references therein): 




Rmdisk V * 'In 5 (3) 



where Mf^sj^ is the disk accretion rate and 77 is a dimensionless factor of order unity. For a star with 
an aligned dipole field, 0.5 < 77 < 1.0. In addition, the mass transfer rate throughout the disk is 
assumed to be comparable to an observational inferred accretion rate from protostellar disks onto 
CTTSs, which decays with time t over the timescale ^ 10^ years and can be fitted by 



—3/2 

where is the Earth mass (Ida & Lin 2004). 

The equation for the spin evolution of the protostar is then given by 

- — - Tdisk)' (5) 

In the above equation, is the moment inertia of the CCTS and T^isk — T^ag + 2acc is the net 
torque arising from the presence of the disk, where Tacc is the torque spinning up the star due to 
the disk accretion: 



Tacc — ^disk-Rmdisk^disk^ 



(6) 



with ^disk being the Keplerian angular speed at the radius Rmdisk (lArmitage fc Clarkd Il996l : 
Matt fc Pudrita l2005l). i^^ and are obtained from the interior structure of the contractin; 



CTTS modeled by lArmitage fc Clarkd (jl996l ). who modified the Eggleton code (jPols et al. 1 11^ 
for the pre-main sequence evolution. 

When the star rotates faster than the inner edge of the disk (i.e. > ftdisk)i the mass 
accreting onto the star should be ceased by the centrifugal force and Tacc is therefore set to be zero 
in our model. Meanwhile, the star experiences immediate magnetic braking as a result of its strong 
fields, leading to disk locking of the stellar spin. Disk locking makes the star and the inner disk edge 
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corotate again and resumes the accretion. Therefore precisely speaking, the disk density near the 
inner edge of the disk, Rmdiski and M^isk should fluctuate with time. In practice, we have simply 
employed equations ([3]) and (jH) to implement the disk properties rather than modelling the detailed 
disk structure near Rmdisk- This is certainly a limitation of our model. With regard to the overall 
disk accretion rate, as we shall see in §3, the fluctuation of the balance between the mass accretion 
and the magnetic braking (i.e. disk locking) is rather fast compared to the timescales of interest, 
meaning that equation (j4j) is a proper approximation to describe the average mass accretion rate 
in our problem. 



2.2. Star-Planet Magnetic Interaction 

The torque on the planet due to magnetic linkage to the star is modelled as the magnetic stress 
{BzB(f)/4:7v) multiplied by the cross-sectional area of the planet's magnetosphere and the distance 
r between the planet and its parent star. It is assumed that the planet has a magnetosphere in 
which the gas behaves like plasma. In addition, the magnetosphere is assumed to be as large as 
the planet's size, so the area of the planet interacting with the stellar field is 2tiR^. Applying the 
star-disk magnetic interaction described in §2.1 to the star-planet case, we obtain the magnetic 
torque on the planet 

B^R^Rp / 1 \ 

Tplanet = ^ \^ ~ ^ ) ' 

where <j) is the orbital angular velocity. The magnetic torque changes with the distance between 
the planet and the central star while the planet moves in an eccentric orbit, so we employ the mean 
torque over one orbit period PorUt- That is, 



fP^rbit BlR^Rl f^orbit / I Q\ 



(8) 



To carry out the above integral, we make the approximation that the orbit does not vary significantly 
during one orbital period; i.e., the orbital angular velocity ^ a^nVl — e^/r^ and r ^ a(l — e^)/(l + 
e cos 0), where a is the semi-major axis, e is the eccentricity, and n is the mean motion of the mutual 
orbit of the planet and the star. Then after integrating and being divided by Porbit^ equation ([8]) 
gives an expression for the mean magnetic torque on the planet: 



d2 d6 d2 
^ -'-planet ^ ^ . 



2 + 3e2 2n, 



a5(l - e2)7/2 a5n(l - e2)2 



(9) 



Here, we have introduced a coefficient e, where < e < 1, to regulate the magnetic torque on the 
planet, a formulism similar to that applied to the magnetic torques between the Galilean satellites 
and Jupiter (Zarka 2007). In this regard, the energy dissipation associated with the magnetic 
torque can be generally considered as the Poynting flux at the magnetopause scaled by e. Hence, 
e = 1 corresponds to the upper limit of the magnetic torque. 
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2.3. Tides on the Planet 



The planet is assumed to have zero obhquity and zero orbital inclination. The gravitational 
tides on the planet affect the evolution of the eccentricity and the spin of the planet. For the sake of 
simplicity, we do not consider the influence of the planet-star magnetic interactions on the spin and 
eccentricity of the planet. Therefore, the rate of the eccentricity change of the orbit is due solely 
to gravitational tides on the planet and the star. The rate contributed from the tidal dissipation 
in the planet alone is given by gp, which can be obtained from (Mardling & Lin 2002) 

5 



81e /M 



*5P 



■p,* 



ll/2(e)0p,^ 
18 



(10) 



where for the tidal dissipation in the star is also shown and will be considered in §2.4. In the above 
equation, /i(e) = (l + 15eV4+15eV8 + 5eV64)/(l-e2)i3/2 and /2(e) = {l + Ze^ + /%) / {l-e^f . 
Moreover, Q'^ and Q'^ are the tidal quality factors of the star and the planet, respectively. All the 
uncertainties associated with the physical processes of tidal heating are contained in the Q' values. 
For a fiducial value, it has been suggested that the Qp values inferred for Jupiter from lo's orbital 
evolution is 5 x 10^ <Qp <2x 10^ (Yoder &Peale 1981). With this Qp value, orbits of planets 
with Mp and Rp comparable to those of Jupiter and with a period less than a week, are circularized 
within the main-sequence lif e span of solar-t ype stars. Q'^ has been inferred to be ^ 1.5 x 10^ for 
the very young binary stars (jLin et al .11 19961 ). 



Since the sy nchronous timescale of the planet is much shorter than the circularization timescale 
(jOu et ahlbood ). we simply assume that the planet's spin equals the orbital rotation throughout the 
evolution. This assumption may be invalid during the phase of the Roche-lobe overflow if the orbit 
changes dramatically due to quick mass loss (see eq. (|T5]) below), causing a large discrepancy be- 
tween the planet spin and orbital motion. This means that the tidal heating due to synchronization 
is unrealistically ignored during this special phase in our model. 



Moreover, the tidal process of orbital circularization shrinks the planet's orbit; i.e., 

2ee 



1 



(11) 



During the orbital circularization, there is substantial tidal heating within the planet. The tidal 
heating of the planet may have contributed sigijnif icantly to the thermal budget that governs the 
planet's physical properties, including its radius flBodenheimer et all I2OOII : IGu et all [20031, l2004l : 
Mardlingl[2007l : llbgui fc Burrowslbood ). 



Ignoring the rotation energy of the planet, we can write the first law of the thermodynamics 
for a tidally heated planet as follows: 



(12) 



where Ende is the tidal heating rate (see Mardling & Lin 2002 for the expression) , C is the intrinsic 
luminosity from the photosphere of the planet, U is the internal energy, and W is the gravitational 
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potential energy. Suppose that U — quGMp/Rp and W = —qwGMp/Rp with the coefficients qu 
and qw depending on the interior structure of the planet. Then in the absence of mass loss (i.e. 
Mp — 0), the above e nerg!;y equation gives rise to the expansion/contraction rate of the planet (cf. 



Mardling fc Lin 



200 



IS 



dRp 



Efide ~ ^ 



M„=0 



(13) 



Even though the expanding/contracting planet is not in the thermal equilibrium, qu, qw, and C can 
be simply expressed in terms of functions of Rp and Mp as long as the timescales under consideration 
are > the convective turnover time for an almost convective planet with the same solid-core mass 
and a thin radiative envelope of the same opacity (i.e. an almost polytropic interior, see Gu et 
al. 2003, 2004). In our approach, the core mass of the planet is assumed to be 19.4M0, which 
is motivated by the core accretion model for the formation of giant planets (Pollack et al. 1996). 
However, the thermal properties of the core are not modelled in our work. qu{Rp), qw{Rp), and 
C{Rp) are obtained only from certain different M p by fitting ; to th e numerical data for the interior 
structures of tidally inflated Jupiters described in lGu et all (|2003l ). Subsequently, their values and 
derivatives with respect to Rp for any desired Mp can be interpolated from these available fitted 
functions. 



Once the planet's radius exceeds its Roche radius Ri (see the Appendix), the mass loss ensues. 
During the planet overflow phase, the Rp adjustment after losing mass plays an important role in 
determining the subsequent evolution. Rp changes with Mp as well according to the equation of 
state of the planet even though the entropy of the planet does not change. Since the interpolated 
data derived from our interior structure code is unable to deal with the processes that are as quick 
as the adiabatic mass loss, we adopt the following simple methodology. Assuming a polytropic 
structure of a planet, the adiabatic mas s-radius e xponent a = {din Rp/ din Mp)s can be estimated 



from the polytropic index of the planet (jGu et al.ll2004 ). a can then be related to the gravitational 
potential energy. It can be easily shown that a = {3 — 4gvi^)/(3 — 2qw)- For simplicity we shall use 
this relation to estimate a from qw for a Roche-lobe filled planet, albeit these polytropic relations 
are derived for an isolated gravitationally bound body. Therefore, the complete expression of Rp 
in the presence of the mass loss through the Lagrangian 1 (LI) point should read 



Rp — 



OR, 
dt 



+ a 



M„=0 



1 1 



(14) 



where we shall adopt the formalism in lKolb fc Ritterl (|l990l ) to estimate Mp (see §2.5 for the detail). 
In the case of a giant planet with large radius such as a young hot Jupiter, a is always negative. 
In addition , \a\ increases w ith i?^, because the planet degeneracy is being lifted when its density 



goes down flGu et al.ll2004l ) 



The internal energy is ignored in iMardling fcTk] ^200i ) 
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Similar to the physical condition for the formation of the magnetospheric cavity of the disk, 
planet's overflow may either form a ring/disk or funnel along the stellar field lines onto the CTTS 
depending on whether Mp is strong enough to dominate over the stellar fields. Therefore, we 
modify equation ^ by replacing M^isk with Mp to define the "overflow cavity radius i^^c" where 
the Reynolds stress of the overflowing gas is comparable to the magnetic stress of the stellar 
fields. Thus, if this radius is larger than planet's semi-major axis, where the overflow occurs, the 
overflowing gas will flow through the magnetic field to the central star. The stellar spin will speed 
up under the action of the torque from planet's overflow, = Mp(GM*a)^/^, which is an analogy 
to the spin- up torque due to the disk accretion described in equation This leads to the evolution 
of the stellar spin caused by the planet's overflow: 

^ = jM,{GM.af'\ (15) 

On the other hand, if the overflow cavity radius is smaller than the semi-major axis, the overflowing 
gas will form a ring/disk around the star instead of directly funnelling along the field lines onto 
the star. In this stage, we assume that the total angular momentum of the overflowing gas is 
immediately transferred to the orbit of the planet, leading to the outward migration of the mass- 
losing planet: 

da\ 2Mpa ^^^^ 



J overflow 

and therefore no spin-up torque is exerted on the star. 

The above equation for outward migration suggests that the mass-losing planet can be detached 
from its Roche lobe as long as Rp does not expand as fast as Ri. Therefore, it is worthwhile to 
examine the linear stability of a Roche-lobe filled planet against mass loss. The definition of Ri 
from the Appendix and equation ([TBI) give Ri/Ri — (— 5/3) (M^/M^). Equation (fHl) then indicates 
that the expansion rate of Rp is ^ —aMp/Mp provided that the thermal adjustment of Rp is less 
important than the adjustment of Rp due to the adiabatic mass loss. Once the planet's photosphere 
reaches the Roche radius (i.e. Rp = i?^), the above analysis infers that the mass loss of a Roche-lobe 
filled planet is linearly unstable (i.e. undergoes exponential growth) if \a\ > 5/3. 



2.4. Tides on the Star 

In addition to the tides on the planet induced by the star, the planet also raises the tides on 
the star to circularize the orbit and to lock the rotation of the star. 

The circularization rate described in the preceding subsection has to be complemented by the 
tidal dissipation in the star; namely, 

^ =9p + 9^, (17) 
where has been displayed in equation (fTUl) . 
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The process of the tidal locking of the CTTS leads to the exchange of the angular momentum 
between the planet's orbit and stellar spin. The stellar spin is assumed to be normal to the orbital 
plane of the planet and t he disk. Owing; to thi s tidal effect, the angular momentum of the stellar 
spin changes at the rate fjMardling fc Linll2002l ) 

5 



J He 



/4(e)f2* - /3(e)n\ / R 



(18) 



_ Mt J \ 2QJ{1 - e2)i/2 

where Jq = Mp[GMta{l — e^)]-^/^ is the total angular momentum of the planet's orbit, /3(e) = 
(l + 15eV2 + 45e''/8 + 5eVl6)(l -e2)-6, and /4(e) = (1 + Se^ + 3e''/8)(l - e^)-^/^. The evolution 
of stellar spin associated with equation (fT8|) is then given by 




[/3(e)n - /4(e) ri* 



(19) 



2.5. Equations for simulation 

After describing the various interactions separately, we list and recapitulate the final equa- 
tions incorporating these interactions for solving the orbital evolution of a young hot Jupiter in a 
magnetospheric cavity. 

The rate of orbital circularization (i.e. eq. (fTTl) ) is given by 

de 

The evolution of planet's size (i.e. eq. (fMll ) has the form 

The first term on the right side is due to thermal adjustment, whereas the second term is due to 
Roche-lobe overflow driven by an adiabatic process that occurs much faster than the first term. 

When Rp > Ri, the mass loss occurs in our simulation. This implicitly implies that we neglect 
the mass loss from the optically thin atmosphere before the planet' s photosphere r eaches its Roche 



radius. We also ignore the possible mass loss through the L2 point (jGu et al.ll2003l ). The mass-loss 
rate of the planet is estimated from the 1-D interior structure of the planet along with the formula 
in Kolb & Ritter 1990 (and see the definitions of the notations therein): 

M, = 2,F,(,)— F3(r:)(-) pd*. (20) 

The orbital decaying rate of the planet is given by 

da _ 2aee 2 < Tpignet > 2aj^ 2Mpa 

dt~ M^anVT^ Jo Mp ' ^ ^ 
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The first term on the right side represents the orbital decay due to circularization (i.e. eq. (fTTI) ). 
The second term is due to the magnetic torque on the planet (i.e. eq. (|9j)) and the third term arises 
from tides on the star (i.e. eq. (jl8l) ). The last term is included only when the overflow occurs (i.e. 
eq. (fTBl) ). Moreover, when Rmdisk shrinks such that the 2:1 reso nance location lies in the disk, 



the planet is expected to slowly migrate inwards mice et al. 1 120081). However, when this happens. 



we assume just for the sake of simplicity that a is decreased instantaneously to maintain the 2:1 
resonance with the inner edge of the disk. 

The stellar spin evolution is determined by a number of factors: 

^ = - Td^sk- < Tpianet > +Mp[GM,a{l - 6^)]^/^} + UJ,. (22) 

The flrst term in the curly bracket arises from the evolution of the stellar moment of inertia during 
the pre-main-sequence phase and the second term represents the torque exerted by the disk on the 
star due to disk accretion and the magnetic linkage (i.e. eq. (ED). The third and the fourth terms 
are due to the magnetic linkages with the planet (i.e. eq. ([H])) and the overflowing gas from the 
overfilled Roche lobe (i.e. eq. (fTSl) ), respectively. The last term cu^ results from the tides on the 
star (i.e. eq. (fTOll ). 



3. Numerical Results 

To solve the ODEs in the preceding section (eq. ([Hj), ([Mj), ([20j), ([ZI]), ([22])), we employ 
the 4th order Runge-Kutta method to study the inward migration of a young hot Jupiter in the 
magnetospheric cavity due to tidal and magnetic interactions with its CTTS. The simulations start 
from t = 10^ yr to 10^ yr. The planet is not introduced to the orbit with a = a2:i until t — 7 x 10^ 
years , a timescale comparabl e to the Type II migration time of a giant planet in a proto-planetary 



disk (jLin fc Papaloizoulll985l ). 



In our model, Sq, e, and Mp are parameterized to investigate the various evolutions of the 
planetary migration in the magnetospheric cavity. Besides, Qj^ Qp and 77 are fixed as 3 x 10^, 10^, 
and 0.75 throughout the simulation. 

The initial period of the central star is set to be 15 days at t = 10^ years (see the left panel 
of Figure Ej). The figure shows the stellar spin will be locked by the disk in a few xlO^ years and 
therefore the initial spin is insensitive to the tidal and magnetic interactions, which proceed on 
the timescales of ^ 10^~^ years. The final disk-lock period (i.e. the Kepleria n period of the disk 



inner edge) depends on the magnitude of the magnetic field of the central star (jArmitage fc Clarke 



19961 ). We discuss two cases of cavity sizes corresponding to Bq = 1500G (§3.1) and = 500G 



(§3.2) separately. The disk with a big cavity will lock the centr al star at the rota t ional period of 



^ 7 days and that with a small cavity will lock it at ^ 3 days. iHerbst fc Mundtl (120051 ) revealed 
the peculiar bimodal distribution of periods for young stars in the Orion Nebula Cluster. The 
distribution of periods at an age of about 1 Myr exhibit two peaks, one located at a period of about 
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8 days and the other at a period of about 2-3 days (jBouvierl 120071 ) . The disk-lock periods in our 
model approximately match the observation. 

We consider the planets of four different planetary masses 2Mj, 1.5Mj, IMj, and 0.7Mj to 
study their different evolutions of migration. From the observations, the close-in transiting giant 
planets have radii ranging from ^0.8 Rj to r^URj. Thus 1.7Rj could be considered as a lower 
limit of the radius of a young Jup iter although some of the large plan et's sizes have been postulated 
to be attributed to tidal heating ( Miller et al. 20091 ) or evaporation (jsaraffe et al.|[20o3 ). Also, 1.7 



Rj corresponds to the theoretical radius of a young Jupiter (Marley et al. 2007). Thus we fix the 
initial Rp = 1.7 Rj for all the cases. 

When t = 7 X 10^ yr and the planet is introduced to the initial location at a = a2:i, the star 
has been locked by the disk. The right panel of Figure [3] shows the decrease in the periods of stellar 
spin and planet's orbit before t < 10^ years. This is caused primarily by the contraction of the 
cavity as i?* decreases (see eq. (|3j)). Orbiting interior to the co-rotation radius, the planet then 
migrates inwards from the 2:1 resonance due to faster revolution of the planet than the spin of 
the central star. This is illustrated by the black curve during the time t > 10^ years in the right 
panel of Figure [3] where the planet with IMj and = 0.16 is used as an example. Meanwhile, 
the star should spin up via the tidal and magnetic interactions with the planet. However, the disk 
locking dominates the evolution of the stellar spin as indicated by the overlap between the red solid 
and blue dashed curves in the right panel of Figure [3l In other words, the inner edge of the disk 
acts as the sink of total angular momentum of the system in our model. In fact, the rotational 
period of the disk inner edge fluctuates slightly from the spin period of the star on short timescales. 
Nevertheless, these two periods are almost the same over long timescales, justifying the argument 
given in §2 that equation (j4j) is a proper approximation to describe the average mass accretion rate 
and determine the cavity size in our problem. 

As mentioned in §2.3, the planet will lose its mass if the planet's radius exceeds the Roche 
radius. However, if the planet's mass decreases to 20 earth masses, the simulation will be forced to 
stop. This is because the core mass of a giant planet is 19.4 earth masses in our model. 



3.1. The Large-Cavity Case 

A large cavity is opened up by the strong stellar magnetic fields with the scaling factor Bq = 
1500 G. In this case, the planet starts to migrate from ^ 0.045 AU, which corresponds to a2:i 
at t = 7 X 10^ years. 

The migration evolutions of the planets with same masses will differ as a result of different ini- 
tial eccentricities. The distance, over which the planet can migrate inward from the initial position, 
is shortest when = and increases as increases. This is because firstly, the circularization of 
a eccentric orbit makes the orbit shrink. Secondly, the inflation of a planet due to tidal heating 
enhances the star-planet magnetic interaction. Also, the tidal torque on the star caused by the 
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planet is more significant as is larger. All of the above factors accelerate the migration rate of 
the planet. When > the critical eccentricity Cc, the planet's radius will exceed its Roche lobe 
and start to lose mass onto the central star. 

We perform the calculations for three cases of ef. (1) = 0; (2) = Cc] (3) = ec — 0.01 for 
each planet mass. The orbital evolutions are shown in Figure [4] and the explanations are presented 
below. 

In the case of = (solid lines in the figure), the results for different masses are not signifi- 
cantly difference. Without an eccentric orbit to bring the planet closer to the central star (i.e. the 
perihelion), the tidal and magnetic effects are excited less effectively in a large distance between the 
planet and the central star. Hence, planets with different masses initially in circular orbits migrate 
to almost the same final positions, a ^ 0.037 AU at t = lO'' yr, as the cavities contract and drive 
the planets to the 2:1 resonances. 

When Ci = ec (dotted lines in the figure), the planets start to overflow their LI points and lose 
mass. As shown in the figure, their evolutions are terminated before the end of the simulations 
{t = 10^ years), suggesting that all of the planets are disrupted. The termination of the simulations 
occurs when the discrepancy between Rp and Ri is so large during the overflow phase that Mp is 
too large to be resolved in the simulation with the time step = 100 years. This time step of 100 
years is already shorter than the time step used for the interpolated data derived from our interior 
structure code. When the runaway solutions happen, \a\ is always larger than 5/3, which is in 
agreement with the linear instability of the mass loss. For the 0.7 and IMj cases, the mass loss 
rate is high. The overflow gas is able to form a ring/disk and hence the planet migrates outward 
according to equation (fT6l) in our model. However, once their overflows occur, \a\ > 5/3 at the 
same time owing to their large Rp. The radius adjustment is predominated by the adiabatic mass 
loss instead of tidal inflation, leading to the runaway solutions to the ODEs. 

On the other hand for the 1.5Mj and 2Mj cases, they migrate fast to the locations closer 
to the central stars owing to the stronger tides raised by these massive planets on their central 
stars. In addition, these more massive planets are more difficult to inflate by tidal he ating due 



to their higher gravitational binding energies and cooling rates (iGu et al.l l2003l . 120041 ) . Hence, 
when they overflow their Roche lobes, Rp is not too large and therefore \a\ < 5/3, leading to the 
stable overflow. Although the overflowing gas initially funnels to the central star, the overflow 
rate becomes substantially high shortly such that Rqc < cl- Consequently the massive planets have 
undergone outward migration almost since the overflow began. The outward migration is primarily 
driven by inward migration due to the tidal and magnetic interactions rather than tidal inflation. 
The overflow finally turns to be unstable when \a\ increases to 5/3 as a and therefore Rp increases. 
The evolution then stops and the planet loses all its gas to the central star. During the phase of the 
fast outward migration, the planet's spin may become unlocked with its orbital motion. Despite 
not being included in our model, the tidal heating associated with the asynchronism would speed 
up the planet inflation and lead to the same catastrophic fate of the planet. Comparing these cases 
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for different initial M^, massive planets require higher Cc thus more tidal heating to be inflated 
against their greater gravitational binding energies and cooling rates. Mp — 2Mj requires ec — 0.35 
and destructs at a ^ 0.02 AU, while Mp = 0.7Mj needs Cc = 0.22 and destructs at a ~ 0.035 AU. 

Now we turn to the case for = ec — 0.01 (i.e. is slightly smaller than Cc) to study the 
conditions in which the planet of a given mass can migrate inward the farthest and still survive in 
the end of the simulation. As can be expected from the previous results for = ec, the dashed 
curves in the figure show that a more massive planet undergoes faster migration due to the stronger 
tidal torque. Thus, the planet with Mp = 2Mj and e^ = 0.34 can arrive at a ^ 0.024 AU, while 
the planet of Mp = 0.7 and 1 Mj can not migrate to the region < 0.03 AU. Furthermore, a low 
mass planet is more easily inflated by tidal heating, giving rise to the relatively stronger magnetic 
interaction (see eq. (Ej)). This explains why the planet of 0.7Mj migrates slightly farther in than 
the planet of IMj in the end of the simulation. 

In equaiton ([9j), the coefficient e regulates the magnetic torque on the planet, and the upper 
limit e = 1 has been applied to obtain the above results. Here, we employ e = to study the orbital 
evolutions in the absence of the magnetic torque. In comparison to the previous orbital evolutions 
shown in Figure [H Figure [5] shows the lower migration rates because of the lack of the star-planet 
magnetic interaction for the Mp =0.7, 1, and 1.5Mj cases. However, the stellar tides induced by 
the 2Mj planet are far dominant over the magnetic interaction, so the massive planet still migrates 
to a ~ 0.016 AU. That position is even farther inward than a ~ 0.024 AU for the e = 1 case as the 
planet with even higher e^ = 0.38 can still survive in the absence of the magnetic torque. Moreover, 
Figure [5] confirms the previous subtle result about the dependence of the magnetic torque on Rp'. 
the slightly more inward migration of the planet of 0.7 Mj than IMj does not reappear because 
these exists no magnetic torque to speed up inward migration for a lighter planet of a larger inflated 
radius. As to the cases for e^ = ec, the simulated planets overflow their Roche radii and perish 
in the same way as in the e = 1 cases; namely, \a\ of various cases finally is > 5/3 and the mass 
loss becomes runaway. For the cases of 0.7Mj, IMj, and 1.5Mj, \a\ > 5/3 once their overflow 
occurs. On the other hand, after reaching the Roche lobe, the planet of 2Mj migrates outward for 
a period of time due to stable overflow and inward migration driven by the tidal torque. It is finally 
destroyed due to the large \a\ that increases with Rp as the mass-losing planet moves outwards. 



3.2. The Small-Cavity Case 

The weaker stellar fields due to the smaller scaling factor Bq = 500 G form a smaller mag- 
netospheric cavity, so the initial position of the planet is much closer to the central star; namely, 
ai 0.032 AU. The planet and the star exert remarkable tidal forces on each other in such short 
distances, and the magnetic torque on the planet is not negligible as well. Ri is getting smaller while 
the planet migrates inwards. As soon as Ri < Rp, the overflow occurs, leading to the destruction. 
Consequently, even e^ = (i.e., no tidal heating in the planet), all the planets destruct irrespective 
of their mass. If Mp is bigger, the migration is faster due to the stronger tidal torque. Therefore, 
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the more massive the planet is, the quicker the planet destructs, as shown in Figure [6] for e = 1. 

The orbital evolutions of the the IMj and 2Mj cases are characterized by the alternating 
out/inward orbital migrations. Once the overflow occurs, the large expansion of Ri due to the high 
mass loss rate dominates over the adiabatic inflation, hence leading to the fast outward migration. 
As a result, Rp temporarily detaches from Ri for a while until the planet moves in and the next 
Roche overflow occurs again. Finally, the planet will lose all its mass when the overflow occurs 
accompanied with high \a\ (i.e. > 5/3). 

In the absence of the magnetic torque (i.e. e = 0), the tides raised by massive planets on the 
stars still drive significant inward migration of the planets. Figure [6] shows that shutting down the 
magnetic torque cannot change the outcome of planet destruction for Mp > IMj. However, the 
planet of 0.7 Mj can survive in 10^ years because the inward migration driven by the weak tidal 
torque from the low mass planet is slow. Therefore, we estimate ec = 0.18 for the planet to perish 
and show the evolution in the leftmost panel of Figure El 

4. Summary and Discussion 

We construct a simple model to study the orbital evolution of a young hot Jupiter in an 
eccentric orbit inside a magnetospheric cavity of a proto-planetary disk around a CTTS. For the 
sake of simplicity, we assume that the magnetospheric cavity of the protoplanetary disk is truncated 
by a stellar dipole field. Through the magnetic linkage between the central star and the disk, the 
rotational period of the star is quickly locked by the inner edge of disk. We then introduce a 
young hot Jupiter at a2:i and restrict ourselves to the orbital evolution of a hot Jupiter inside 
the cavity. Assuming that the planet behaves like a plasma, we apply the same formulism for the 
star-disk magnetic interaction to that for the star-planet interactions. We adopt the equilibrium- 
tide equations with = 3 x 10^ and Qp = 10^ to model the star-planet tidal interactions. In so 
doing, our model focuses on the planet migration due to tidal and magnetic interactions between 
the planet and the star, and does not consider any interactions between the planet and the disk. 
However, as the size of the magnetospheric cavity evolves, the planet is artificially pushed to a2:i if 
it lies beyond a2:i. We vary three parameters Bq (= 500 and 1500 G), initial Mp (= 0.7, 1, 1.5, 2 
Mj), and initial e (e^ = 0, — 0.01, and Cc) in our simulations to investigate the fate of the planet 
migration under the influence of the cavity size, planet mass, and orbital eccentricity. 

Two sizes of magnetospheric cavities are considered to approximately match the bimodal dis- 
tribution of spin periods of young stars in the Orion Nebula Cluster. The initial a2:i corresponding 
to the large (small) cavity opened up by the stronger (weaker) stellar fields is ^ 0.045 AU (~ 0.032 
AU). In the case of the large cavity {Bq = 1500 G), the planets require nonzero to enhance their 
tidal and magnetic interactions with their central stars and make significant inward migration. The 
migration rate increases with for a given planet mass. When is as large as the critical value 
Rp exceeds Ri and the Roche overflow occurs. For planets with different masses, massive planets 
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{Mp > IMj) have higher ec (i.e. more intense tidal heating) than less massive ones {Mp < Mj) 
because massive planets are more difficult to thermally inflate due to their greater gravitational 
binding energies and cooling rates. As a result, stronger stellar tides raised by the planet with 
higher Mp and allow for faster migration. Overall, massive planets can therefore migrate further 
in than less massive ones, as shown in the simulations for = — 0.01. When = ec, low-mass 
planets {Mp = 0.7 and 1 Mj) inflate by tidal heating faster than inward migration. Hence they 
overflow their large Roche lobes at large a, resulting in the low density and therefore the runaway 
mass loss. On the other hand, the high-mass planets {Mp = 1.5 and 2 Mj) migrate inwards fast 
without significant tidal inflation. In contrast to the low-mass planets, once these planets fill their 
small Roche lobe at small a, their density is high enough that they can undergo stable LI overflow. 
During this phase, tidal and magnetic interactions instead of tidal inflation drive them to lose mass, 
expand, and migrate outward. This stable mass-loss phase proceeds until these outward migrating 
planets are large enough to become unstable against mass loss. Finally, they lose all gas to their 
central stars. 

In the case of the small cavity {Bq = 500 G), the simulated planets in circular orbits all quickly 
migrate in due to the fierce tidal and magnetic interactions until they stably overflow their Roche 
radii. After that, they can move outwards due to the mass loss and move inwards again by the 
tidal and magnetic torques. As their mass goes down and their degeneracy is lifted, all of these 
planets are finally destroyed, suffering from the runaway mass loss as in the large-cavity case. 

To study the significance of the magnetic interactions between the star and the planet in our 
model, we also simulate the cases by shutting down the interaction (i.e. e = 0). We found that the 
migration of less massive planets is more sensitive to the magnetic interaction, which is enhanced 
by their easily inflated radii. Setting e = makes the migration rate of low-mass planets slower, 
as shown in the case for Mp = 0.7Mj with = 0.21 in the large-cavity case. By contrast, the 
inward migration of massive planets is less sensitive to the magnetic interaction but more to the 
tidal interaction. As has been summarized in the preceding paragraph, all planets in the small 
cavity destruct under both the tidal and magnetic torques at such close-in initial positions from 
the stars even if = 0. When e = 0, only the planet of 0.7 Mj can survive unless the low-mass 
planet starts with an eccentric orbit with ec = 0.18. 

The mo del presented he re for the orbital e v olutio n of young hot Jupiters is different from 
the works by iGu et al.l (120031 ) and iTrilling et al.l (119981 ) . In this work, we empl oy simple mode ls 
including the star-planet interactions and disk locking, which are not considered by lGu et al.l (120031 ) . 
Furthermore, we take into account the pla net expansion due to the adiaba tic mass l o ss and model 
the mass loss rate using the approach in iKolb fc Ritterl (Il99d ) , whereas IGu et al.l (120031 ) ignore 
the radius adjustment to the mass loss and use a free parameter to control th e mass lo s s rate . 
Ignoring the mass-radius relation during the Roche-lobe overflow phase as done in lGu et al.l (120031 ). 
a Roche-lobe filled giant planet with = Cc continuously loses its mass via tidal inflation and 
migrates outwards until the tidal heating rate is weaker than the cooling rate, leading to a survived 
planet of lower mass. On the other hand. ITrilling et al.l (119981 ) consider young hot Jupiters moving 
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inwards via the Type II migration without a magnetospheric cavity and tidal inflation. In their 
work including the mass-radius relation, once the planets of mass < 3.36Mj get close enough to 
their central stars and therefore overflow their Roche lobes, they migrate outwards and finally lose 
all their mass onto the star. In this regard, our results are compatible with theirs. On e of the major 
differe nces is that the mass loss is driven by the Type II migration in the model by i Trilling et all 
(119981 ). while in our work the mass loss is driven by the loss of orbital angular momentum via 
the tidal and magnetic interactions with the parent star whose spin is almost locked by the disk. 
Moreover, the planets of mass < 2Mj with ei < Cc survive in the cavity during the CTTS phase in 
our model. 

For the last few years, transit surveys have revealed the peculiar Mp — a correlation for hot 
Jupiters: less massive hot Jupiters (< 1 Jupiter mass) are almost absent within ^ 0.03 AU from 
their parent dwarf stars of the ages ^ Gyrs. Our results show that during the CTTS phase, planets 
of Mp > IMj with modest initial eccentricities (e^ ^ 0.3) and sufficient magnetic interactions in 
a large cavity can migrate to a < 0.03 AU, while planets of Mp < IMj cannot. More specifically, 
given ei < ec, the planet of 2Mj can safely arrive at a < 0.024 AU and the planet of 1.5 Mj can get 
to a < 0.025 AU in the end of our simulations. Whether or not these planets can further migrate 
toward their spun-down parent stars during the main-sequence phase depends on Q'*, which is not 
yet well understood. 



Ogilvie fc LinI (120071 ) suggested that induced by hot Jupiters may become large (i.e. » 10^ 



as CTTSs evolve to ma i n-sequence stars based on their dynamical-tide model. As they and others 
(jPatzold fc Rauerll2002l : I Jiang et al.ll2003l ) pointed out, we would be extremely fortunate to observe 
these very close-in hot Jupiters if the current Q'^ ^ 10^, implying that of the solar- type main- 
sequence host stars may be quite large. However, it should be noted that these very close-in planets 
were all discovered by transit surveys (se e Figure [1]) . Although transit surveys overall prefer to 
detect close-in planets (iGaudi et al.ll2005l ). the detectability of a hot Jupiter as a function of the 
semi-major axis of interest here (i.e. a < 0.05 AU) is extremely ill-defined as the various transit 
surveys are subject to different obse rvational limitati ons such as the transit depth, sampling rate, 
red noise, and detection threshold (jPont et al.ll2006l ). If indeed gradually becomes » 10^ as 
CTTSs evolve to main-sequence stars, therefore allowing the planets to move a little bit inward 
during the short transition after t = 10^ years while halting the further planet migration for most 
of the main-sequence phase, our model may have the potential to explain the absence of low-mass 
hot Jupiters within a < 0.3 AU. 

On the other hand, a large body of studies have suggested that Q'^ associated with the tidal 
dissipation induced by the hot Jupiters in main-sequence dwarfs could be as small as ^ 10^~^. As a 
result of the efficient tidal dissipation, the orbits of the hot Jupiters can decay so significantly that 
the planets niay pl unge into their host stars during the main-sequence phase (jjiang et al.l l2003l : 
Jackson et al.ll2008l ). This "accretion" model may acc ount for the pauc ity of the extremely close-in 
hot Jupiters concluded from radial- velocity surveys (IGaudi et al.ll2005l . also see Figure [1]). In addi- 
tion, the relatively young ages of the extremely close-in hot Jupiters revealed by transit surveys may 
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also lend support to the model (I Jackson et al.ll2009l ). Although the planet-met allicity correlation 



indep endent of the spectral type does not seem to favor the accretion model (IFischer fc Valenti 



20051 ) , excessive heavy elements in the shallow convection zone of high-mass dwarfs may be able to 
sink down to the radiative zone via the double diffusiv e effect, thereby eliminating th e metallicity 
enhancement and resolving the problem for the model (I Vauclairl |2004|) . Furthermore, |Pont| (120091 ) 
found an excess spin of the host stars of hot Jupiters, which may arise from tidal spin-up, although 
the uncertainty of the stellar ages and the limitation of a small number of sam p les de mand more 
detections in the future to confirm the conclusion. The study by Ijackson et all (120091 ) for the ac- 



cretion model is carried out based only on one-mass case (i.e. IMj). If the accretion model is the 
main cause responsible for the observed Mp — a relation, it would be expected that the extremely 
close-in giant planets of mass < IMj of younger ages should also be detected by transit surveys. 
The orbital decay timescale during the main-sequence phase due only to the tidal dissipation is 
given by tude ^ (-^o)/(2J*) (see equation (fT8]) ). In the case of a planet of 0.7Mj, the planet can 
survive in ^ 10^ yrs from the distance of ^ 0.03AU at t = 10 Myrs. In contrast, the orbit of a 
planet of 2Mj will decay on the timescale of ~ 10^ yrs from a distance of ~ 0.025AU at t = 10 
Myrs. In other words, during the main-sequence phase, less massive planets excite weaker tides 
on their parent stars and therefore still stay outs ide a ^ 0.03 AU, while the n iassive planets kee p 
falling into the stars in the way as interpreted bv lPatzold fc Rauerl (120021 ) and I Jiang et al.l (120031 ). 
The recently discovered hot Jupiter WASP-18b with an orbital period of 0.94 days a nd a mass of 
lOMj may provide a contraint on Q^^ of a main-sequence star in the next few years (jllellier et al.l 
200j. 



In this paper, we have attempted to bring several physical components to the context of the 
evolution of young hot Jupiters. Nonetheless, a large number of assumptions have been made to 
simplify the physical processes included in our model. The structure of the truncated disk is not 
modelled. How a planet mig rates to a2:i and how the eccentricity evolves during its entry to the 
cavity (e.g. IRice et al.l l2008l ) are not addressed. The disk accretion rate is a ssumed to be quas i- 
steady to determine the cavity sizes, while observations show otherwise (e.g. iBaraffe et al.ll2009l ). 
The dipole configuration of the stellar fields connecting to the disk ignores any types of winds even 
though the effect of disk locking is introduced in a simplified manner to imitate the sink of the total 
angular momentum in our model. The star-planet magnetic interactions are parameterized to scale 
with the Poynting flux at the magnetopause of the planet rather than appealing to any spec ific 
models, such as the dissipation due to the Alfven waves and unipolar inductor (e.g. Izarka 2007 ) or 
induced by tilted stellar dipole fields (e.g. Laine et al. 2008; cf. Papaloizou 2007). The values of 
Q'^ and chosen for the calculations are quite arbitrary; we have not explored a wide range of the 
values. The interior structure of an inflated pl anet and the mass-radius relation of a ho t Jupiter 
are inferred from the 1-D numerical simulation (IBodenheimer et al.ll200ll : iGu et al.ll2003l ) even for 
a Roche-lobe filled planet. Also, the code ignores thermal properties of the solid core of a giant 
planet, leaving a question as to what would happen to the core under the intense tidal h eating. In 
a hot Jupiter-star system, the mass loss may occur via the Lagrangian 2 point as well (jOu et al.l 
20031 ) . which is not considered in the present work. Besides, the mass loss from a planet in an 
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eccentric orbit may be intermittent and nonconservative, while the mass loss rate in our calculation 
is simply estimated from the LI overflow at r = a. To use a more appropriate expression for 
the Roche radius and to consider the m ore realistic mass loss in asynchronous eccentric planetary 
systems, a more careful treatment (e.g. ISepinsky et al.l 120071 ) is desired to refine the results. In 
light of these limitations and possible uncertainties, the simulations covering a wide range of the 
parameter space along with more realistic modelling on individual issues will be explored in a future 
work. 
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A. Appendix 



The Roche radius Ri is defined as the distance between LI and the planet's center of mass. 
The conventional express io n of Rj is derive d from the circular restricted three-body problem (e.g. 
Murray fc Dermottlll999l ). IGu et all (120031 ) assumed that the Roche lobe overflow occurs around 
perihelion and made an approximation that Ri = (M^/3M*)^/^a(l — e). Here we investigate 
the validity of the approximati on by employin g; the expression for the LI point in the elliptical 
restricted three-body problem (ITodoran 1 1 19931) . The exact values of the Roche radius in units 



of a are calculated using the equations in iTodoran I (119931 ) and are shown in Table [T] for various 
mass ratios and eccentricities relevant to our problem. The table shows that the Roche radius 
at perihelion is larger than that at aphelion in some of the cases. It is because the planet orbits 
faster (slower) than the "fictitious" circular motion performed at perihelion (aphelion). As a result, 
the centrifugal force due to the orbital motion is larger (smaller) at perihelion (aphelion) than the 
centrifugal force of the fictitious circular motion at the point, therefore enlarging (reducing) the 
Roche radius relative to the "Roche radius" associated with the circular motion. 

Furthermore, the above analysis is based on the assum ption that the spin of the planet is 
instantaneously synchronized with the orbit at any moment (jPratt fc Strittmatterl Il976l ). While 
it is true for a synchronized circular orbit, in an eccentric orbit the planet's spin is normally 
"synchronized" with i ts orbital motion close to perihelion where the tidal forcing is expected to 
be near its maximum llvanov fc Papaloizoul 120071 : iLangton fc Laughlinll2008l ). Following 

the same line of the above argument, a "synchronized" planet spins faster than its orbit at aphelion, 
implying that the Roche radius is smaller than the value associated with the circular motion at the 
point. 

It should be noted that equation (l20l) for the mass loss rate is also derived for a circular orbit. 
The equation should be modified as the LI point changes its location along an eccentric orbit. 
Besides, when the planet's photosphere is close to its time- varying Roche lobe in an eccentric orbit, 
the outer part of the planet should expand or contract with the Roche lobe on the local dynamical 
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timescale, which introduces additional complexity to estimate the mass loss rate. 

It is for the above complex reasons that we simply use the expression Ri = (Mp/3M*)V3 

a to 

specify the Roche radius of a planet in an eccentric orbit in the present work. 
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Fig. 1. — The semi-major axis vs. mass relation for the hot Jupiters of mass between 30 Earth 
masses and 2 Jupiter masses. The data are adopted from the website |http://exoplanet.eu/[ The 
transiting exoplanets are shown by the asterisks and the exoplanets discovered in radial velocity 
surveys are denoted by the black circles. The exoplanets that are known to be in retrograde orbits 
are excluded in the plot because their orbital orie ntations are unlikely to be explained by the 



planet migration model as presented in this paper (INagasawa et ahl 12008 



Narita et al 



Anderson et al. 



20091; 



20091). 




Fig. 2. — A schematic illustration of our simple model. 



-26- 




Fig. 3. — The evolution of the stehar spin (red hne), the rotational period of the inner edge of the 
disk (blue dashed hne), and the orbital period of the planet (black line). The left panel displays 
that the star's initial spin period (i.e., 15 days) at t — 10^ years is locked to 8 days by the disk 
in ^ 4 X 10^ years. Afterwards, the star corotates with the disk inner edge, and its spin rate is 
determined by the disk accretion and the disk locking. The planet's orbital period is added in the 
right panel to compare to the stellar spin period and the rotational period of the inner edge after 
the planet is introduced ait — 7x 10^ years. The difference between the planet's orbital and stellar 
spin rates leads to the planet's inward migration as shown by the black curve, while the star does 
not spin up via the tidal effect and star-planet magnetic linkage due to the strong disk locking. 
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Fig. 4. — Orbital evolutions of the planets in the large cavities. By modulating the initial eccentric- 
ity, the orbital evolution will change. Each panel displays the results for a different initial Mp. The 
solid line demonstrates the orbital evolution of the planet with = 0; the dashed line illustrates 
the case in which the planet of a given initial mass can migrate the longest distance; the dotted 
line shows that the planet with = Cc destructs at the breaking point. 
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Fig. 5. — The same as Figure H] except that the coefficient e is reduced to 0, meaning that the 
magnetic fields do not exert any torque on the planet. The critical eccentricities Cc are relatively 
high compared to the previous cases for e = 1. 
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Fig. 6. — Comparision between the cases with the magnetic torque in fuh strength (e = 1 plotted in 
sold line) and without the magnetic torque (6 = plotted in dashed line) for the orbital evolutions 
of the planets of differernt masses with = in the small cavities. Note that since the planet of 
0.7Mj with ei = survives up to the end of the simulation at t = 10'' yr, its ec is estimated and 
the corresponding evolution is also shown (dotted line) in the leftmost panel. 



